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The inelastic neutral reaction of neutrino on ''He is calculated microscopically, including full final 
state interaction among the four nucleons. The calculation is performed using the Lorentz integral 
transform (LIT) method and the hyperspherical-harmonic effective interaction approach (EIHIf), 
with a realistic nucleon-nucleon interaction. A detailed energy dependent calculation is given in 
the impulse approximation. With respect to previous calculations, this work predicts an increased 
reaction cross-section by 10% — 30% for neutrino temperature up to 15 MeV. 
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The interest in neutrino reactions with nuclear targets 
stems from the role they play in major questions of con- 
temporary physics. Such reactions are of central impor- 
tance in various astrophysical phenomena, such as super- 
nova explosion and the nucleosynthesis of the elements. 
In this letter, we present a microscopic ab-initio calcula- 
tion of the neutral inelastic reactions of ^He with Vxii^) 
{x = e,fj.,T). 

Core collapse supernovae are widely accepted to be a 
neutrino driven explosion of a massive star. When the 
iron core of a massive star becomes gravitationally unsta- 
ble it collapses until short-range nuclear forces halt the 
collapse and drive an outgoing shock through the outer 
layers of the core and the inner envelope. However, the 
shock loses energy through dissociation of iron nuclei and 
neutrino radiation, and gradually stalls, it becomes an ac- 
cretion shock. It is believed, but to date not proven, that 
the shock is then revived as neutrinos emitted from the 
collapsed core (the proto-neutron star) deposit energy in 
the collapsing layers to overcome the infall and eventu- 
ally reverse the flow to an outgoing shock which explodes 
the star. Hydrodynamic simulations of a collapsing star, 
which are restricted to spherical symmetry, fail in reviv- 
ing the shock Lately it was shown that even in 
full 2-D calculations the shock is not revived. In order 
to revive the shock, the neutrinos must deposit about 
1% of their energy in the matter behind the shock. The 
latter, which is assumed to be in thermodynamic equi- 
librium, is composed mainly of protons, neutrons, elec- 
trons, and '*He nuclei. In contrast to the fairly known 
cross-sections of neutrinos with electrons and nucleons, 
the interaction of neutrinos with "^He is not accurately 
known, and to date there is no realistic microscopic cal- 
culation of the inelastic ^He-neutrino cross-section. The 
effect of neutrino-^He interaction on the delayed shock 
mechanism was investigated by Bruenn and Haxton 
through a presupernova 1-D model of a 15 Mq star. In 
that model, they found only a small reheating of the mat- 
ter behind the shock, which can be attributed to the low 
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mean energy of the neutrinos in comparison to the high 
threshold energy of the Alpha nucleus. This conclusion 
may change with different progenitor, or with enlarged 
inelastic neutrino-^He cross-sections. 

The neutrinos migrating out of the proto-neutron star 
are in flavor equilibrium for most of their migration. 
The electron-neutrinos remain in equilibrium with mat- 
ter for a longer period than their heavy-flavor counter- 
parts, due to the larger cross sections for scattering of 
electrons and because of charge current reactions. Thus 
the heavy-flavor neutrinos decouple from deeper within 
the star, where temperatures are higher. Typical calcu- 
lations yield temperatures of ~ lOMeV for /i- and r- neu- 
trinos jj] , which is approximately twice the temperature 
of electron-neutrinos. Consequently, there is a consider- 
able amount of f^_r with energies above 20 MeV that can 
dissociate the ''He through neutral reaction. 

The flux of neutrinos emitted in the collapse process is 
sufficiently large to initiate nucleosynthesis in the over- 
laying shells of heavy elements. Neutral reactions of Al- 
pha and neutrino in the inner Helium shell are part of 
reaction sequences leading to the production of the rare 
A = 7 Lithium and Beryllium isotopes [3, 0. Thus, 
better understanding of the v — a reaction can lead to 
better prediction for the abundances of these elements. 

Theoretical understanding of neutrino-nucleus scatter- 
ing process is achieved through perturbation theory of 
the weak interaction model. The nuclear electroweak 
transition operator consists of one- and many-body com- 
ponents. The many-body currents are a result of meson 
exchange between the nucleons, and usually contribute 
up to 10% of the cross-section, in the supernova energy 
regime. However, when leading one-body terms are sup- 
pressed their contribution can be even larger. The cur- 
rent work is done in the impulse approximation, thus 
taking into account only one-body terms. The one-body 
currents connect the *He ground state and final state 
wave functions. In order to calculate the cross-section in 
a percentage level accuracy, one needs a solid estimate 
of these wave functions. Alas, for nuclear systems with 
more than three constituents, where particle correlation 
plays a decisive role, the computation of intermediate- 
energy continuum wave function is currently out of reach. 
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TABLE I: Flavor and temperature averaged inclusive in- 
elastic cross-section and energy transfer cross-section cal- 
culated. The temperatures are given in MeV, the cross- 
sections in 10"*'^ cm'^, and the energy transfer cross-sections 
in lQ-^°cm^MeV 



To facilitate the calculation of the neutral reaction of 
neutrino and alpha particle we introduce several modern 
methods. The calculation of the nuclear dynamics is car- 
ried out by combining two powerful tools: the Lorentz 
integral transform (LIT) method and the effective in- 
teraction hyperspherical harmonics (EIHH) method 
First we use the LIT method in order to convert the scat- 
tering problem into a bound state like problem, and then 
the EIHII method is used to solve the resulting equa- 
tions. Using this procedure we solve the final state inter- 
action problem avoiding continuum wave functions. This 
method was used successfully to calculate the photoab- 
sorption cross sections of up to six body nuclei 0- To 
this end we use nuclear Hamiltonian consists of the real- 
istic nucleon-nucleon potential AV8'. 

In the limit of small momentum transfer (compared to 
the Z particle rest mass), the effective Hamiltonian can 
be written as 



where G is the Fermi weak coupling constant, j^{x) 
is the leptonic current, and is the hadronic current. 
The matrix element of the leptonic current is {f\jfi\i) — 
lfj.e~^^'^ , where = 'S(fc^/)7^,(1 — ^^)u{ky). The nuclear 
current, 



jhadromc ^ _ 3 • sin^ ) ^ ^ - 2 • sin^ Ow^J^., 

(2) 

consists of one body weak currents, but also many body 
corrections due to meson exchange. In this work we 
use the impulse approximation. Since the momentum 
transfer relevant to our calculation are small compared 
to the nucleon mass, we ignore relativistic corrections. 
The differential cross-section is given by Fermi's golden 
rule. Thus, in order to consider recoil effects, and with 
unoriented and unobserved targets, the differential cross- 
section takes the form, 
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f helicities 



where kf is the momentum of the outgoing neutrino, uj 
is the energy transfer, and q is the momentum transfer. 

Choosing the z direction to be parallel to the momen- 
tum transfer, and 9 to be the angle between the incoming 
neutrino direction and outgoin g n eutrino direction, the 
cross-section can be written as [l(l |. 
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the — (+) is for neutrino (anti-neutrino). The functions 
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are the response functions with respect to the transition 
operators Oi and O2 (when Oi — O2 we use the nota- 
tion Rq = Rq^). I ^0,/) and Eqj are the wave function 



and energy of the ground and final state, respectively. 
The transition operators Cj{q), Lj{q), Ej{q), Mj{q) are 
the reduced Coulomb, longitudinal, transverse electric 
and transverse magnetic multipole operators. Since the 
relevant energy regime is up to ~ 60 MeV, the main 
operators contributing to the inelastic cross-section are 
the axial vector operators £^2,^2; -^^1 ; ^0 ^^'^ vector 
Ci, El, Li. Usually, the main contribution comes from 
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FIG. 1: relative error in the sum-rule of the leading response 
functions with respect to the hyper-angular momentum quan- 
tum number K. The error bars reflect the uncertainty in 
inverting the LIT. 



the Gamow- Teller operator but due to the closed shell 
character of the ^He nucleus, it is highly suppressed. In 
this energy range the long wavelength limit jl3 is accu- 
rate up to about 5%, thus it is highly informative to look 
at the the long wavelength expansion of these operators 

C'lAiiq) = Fv^YiMir 

Eui{q) = -V2-CiM{q 

q 

LiM{q) = — CiM{q 
q 



qr _ 



M!Miq) = FA^a-YiiM{r 
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qr 
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Here r is the nucleon's location relative to the system's 
center of mass. The response functions are calculated by 
inverting the Lorentz integral transforms 



duj- 



[uo - anY + erf 
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where a ~ (Jii + and | ^j) {i = 1, 2) are solutions of 
the Schrodinger like equations 

iH-Eo-a)\ ^,{a))=6, \ 

The localized character of the ground state, and the 
imaginary part of a, give these equations an asymptotic 



FIG. 2: Temperature averaged inelastic cross-sections at 
temperature T = 10 MeV. The solid line is the differential 



cross-section, (^)t = 
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dashed line is the differential energy transfer cross-section, 
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boundary condition similar to a bound state. As a result, 
one can solve these equations using the EIHH method. 
In this approach, the potential is replaced by an effective 
potential constructed via the Lee-Suzuky method 
and the new equation is solved by expanding '^q and ^'j 
in four-body anti-svmmetrized Hyperspherical Harmon- 
ics basis functions We calculate the matrix element 
(^i I ^2) using the Lanczos algorithm p^ . 

The combination of the EIHH and LIT methods brings 
to a rapid convergence in the Response functions. In 
Fig. ^ , one can see the relative error in the sum-rule of 
the main response functions with respect to the hyper- 
angular momentum quantum number K. It can be seen 
that upon convergence the relative error is well below 1%. 
The error bars presented reflect the error in inverting 
the LIT. Bearing in mind that the cross-section, up to 
kinematical factors, is the sum of the response functions, 
this is a measure of the accuracy in the calculation of the 
cross-section. 

It is well known that realistic 2-body NN potentials 
lead to an under-binding of about 0.5 — 1 MeV for the '^He 
and the triton nuclei and an under-binding of about 3 — 4 
MeV for ''He. For the AV8' force with a simple Coulomb 
interaction we obtained a binding energy of 25.19 MeV 
for ^He, and 7.76 MeV for the triton. Thus our model has 
a discrepancy, A « 2.4 MeV, with respect to the exper- 
imental inelastic reaction threshold. In order to correct 
for this difference we shifted the response function to the 
true threshold, i.e. R{uj) — > R{uj — A). 

It is assumed that the neutrinos are in thermal equi- 
librium, thus their spectrum can be approximated by the 
Fermi-Dirac distribution with characteristic temperature 
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T. As a result, the interesting quantities are the temper- 
ature averaged cross-section and energy transfer cross- 
section: 



d{auj)T d{a)T 
doj duj 



(7) 
(8) 



where f{T,k) is normaUzed Fermi-Dirac spectrum with 
zero chemical potential, temperature T, and energy k, 
i.e., 



0.5546 fc2 

T3 e*^/^ + 1 ■ 



(9) 



As a typical example we present in Fig.|21the calculated 
cross-section for T = 10 MeV. In Table HI we present 
the calculated total temperature averaged cross-section, 
(ct)t = Wi'^v + o'-)t, and energy transfer cross-section, 
{auj)T = Wi^f^u + ^^o'f)t, as a function of the neutri- 
nos' temperature. Also presented are earlier results by 
Woosley et. al. 0. It can be seen that the current 
work predicts an enhancement of about 10% — 30% in 
the cross-section. 

The energy transfer cross-section was fitted by Haxton 
to the formula ^ , 



(crujjT = a 



^ lOMeV J 



(10) 



with the parameters a = 0.62 • 10 ^°cm^MeV, Tq — 
2.54MeV, (3 = 3.82. A similar fit to our results yields 
a = 0.64 • lO-^^cm^McV, Tq = 2.05MeV, /3 = 4.46. It 



can be seen that the current work predicts a stronger 
temperature dependence of the cross sections. For exam- 
ple, a 15% differnce between these calculations at T = 10 
Mev, grows to a 50% difference at T = 16 MeV. 

In conclusion, a detailed realistic calculation of the 
inelastic neutrino-'*IIe neutral scattering cross-section is 
given. The calculation was done in the impulse approx- 
imation with numerical accuracy of about 1%. The dif- 
ferent approximations used here should result in about 
10% error, mainly due to many body currents, which are 
not considered in the current work. 

The effect of these results on the supernova explosion 
mechanism should be checked through hydrodynamic 
simulations, of various progenitors. Nonetheless, it is 
clear that our results facilitate a stronger neutrino-matter 
coupling in the supernova environment. First, our calcu- 
lations predict an enhanced cross section by 10% — 30% 
with respect to previous estimates. Second, we obtained 
steeper dependence of the energy transfer cross-section 
on the neutrino's temperature. Thus, supporting the ob- 
servation that the core temperature is a critical parame- 
ter in the explosion process. It is important to notice that 
the energy-transfer due to inelastic reactions are 1 — 2 or- 
ders of magnitude larger than the elastic reactions, ergo 
the inelastic cross-section are important to an accurate 
description of the Helium shell temperature. 
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